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A mixed integer linear program is presented for deterministically scheduling departure 
and arrival aircraft at airport runways. This method addresses different schemes of 
managing the departure queuing area by treating it as first-in-first-out queues or as a simple 
parking area where any available aircraft can take-off irrespective of its relative sequence 
with others. In addition, this method explicitly considers separation criteria between 
successive aircraft and also incorporates an optional prioritization scheme using time 
windows. Multiple objectives pertaining to throughput and system delay are used 
independently. Results indicate improvement over a basic first-come-first-serve rule in both 
system delay and throughput. Minimizing system delay results in small deviations from 
optimal throughput, whereas minimizing throughput results in large deviations in system 
delay. Enhancements for computational efficiency are also presented in the form of re- 
formulating certain constraints and defining additional inequalities for better bounds. 


I. Introduction 

A IRCRAFT departing from an airport face numerous constraints in the scheduling of their departure times. 

These constraints include wake vortex separation for successive departures, miles-in-trail separation for aircraft 
bound for the same departure fixes, and time-window or prioritization constraints for individual flights. However, 
departure runway operations also include runway crossings by arrival aircraft; arrival aircraft may need to cross the 
departure runway to come to the terminal. Efficient scheduling of departure operations would require scheduling 
runway crossings and departures simultaneously. Addressing all the above constraints in a single framework is 
critical to increasing airport capacity and efficiency for the Next Generation Air Transport System (NextGen) 
concepts. 

Prior research on runway scheduling has addressed some of the above issues, though most of the work has been 
devoted to scheduling arrival aircraft. The aircraft landing problem has been treated as an example of the generic 
job-shop scheduling problem and a dynamic programming algorithm for single machine scheduling was developed 
and applied to aircraft landing 1 . Other papers also address the aircraft landing problem as a job-shop scheduling 
problem 2 " 4 using a mixed integer formulation. Constraint Position Shifting (CPS) has been proposed for the 
dynamic scheduling of arrival aircraft along with heuristics to solve it 5 . The idea of CPS was further incorporated 
into a dynamic programming approach for scheduling aircraft landings 6 and was later extended to departure 
scheduling 7 . CPS reduces the solution space by limiting the number of positions an aircraft can occupy in the 
sequence and thus could lead to sub-optimal solutions. The dynamic programming approach in Ref. 7 has been 
extended into a generalized dynamic program for solving the departure scheduling problem 8 . This method 
determines pareto-optimal solutions for multiple objectives but does not assign aircraft to available “queues” (when 
the taxiway area next to a departure runway is operated as queues). The decision factors behind such pre-determined 
queues are not explicit and could potentially be based on a subjective controller input. In a recent paper, the authors 
determine runway queues for individual aircraft and schedule departures within the same framework, with the 
framework being generic to both with and without queue situations 9 . However, there is little evidence of a unified 
approach for managing departure queues, departure scheduling and runway crossing scheduling at the same time in 
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the current literature. Runway queue assignment, departure scheduling and runway crossing scheduling might be 
done in multiple steps by the same controller. In such cases, the large number of decisions to be made and the 
emphasis on safety might result in poor schedules due to divergence in the three processes, especially in the case of 
increased traffic. There is no evidence to suggest a uniform manner for these processes in current operations, 
although there is evidence of use of a first-come-first-serve rule being used at other parts of an airport 10 . Thus, 
determining runway queues for individual aircraft and scheduling departures as well as arrival crossings within the 
same, generic framework would potentially improve efficiency and throughput of overall surface operations at busy 
airports. 

This paper presents a Mixed Integer Linear Program (MILP) for incorporating active runway crossings in 
departure scheduling, extending on prior work in Ref. 9. The model assigns aircraft to available departure queues 
and schedules runway departure times as well as crossing times, explicitly considering constraints for wake 
separation and departure fix restrictions. The paper builds on previous work on MILP -based departure scheduling 9 , 
and considers the case where aircraft cross from pre-determined queues based on pilot selection. 

The rest of this paper is organized as fol lows : Sectionflfldiscusses the approach, illustrating it using Dallas-Fort 
Worth International Airport (DFW). Section [Tn|p resents the mathematical formulation, describing the notation and 
the interpretation of the expressions. This is followed by a brief note on the computational enhancements. SectionQv] 
compares the results from the MILP -based solution to a first-come-first-serve scheme in a deterministic setting for 
numerous randomized scenarios based on real-world settings. The choice of objective function is also discussed. 
Section[v~|concludes with directions for future research. 

II. Approach 

The MILP takes as input the sequence of departure (and arrival) aircraft for a runway, along with their earliest take- 
off (and crossing) times and an optional prioritization scheme based on a time-window of take-off (and crossing) for 
each aircraft. The model then assigns departure aircraft to the available departure queues and schedules runway 
usage times for both departure and arrival crossings. 

The MILP approach assumes a deterministic entry time for each aircraft at the entrance of the runway queue 
system or crossing, depending on whether it is an arrival or departure aircraft. The holding area adjacent to the 
departure runway is modeled as a series of queues operated on a first-in- first-out (FIFO) basis. The number of 
queues depends on the airport and runway configuration and is an input to the model. This is illustrated in Figure 1, 
using runway 17R at DFW airport as an example. The figure shows multiple taxiways feeding into three different 
“departure queues” which then feed into the departure runway. 

As shown in Ref. 9, there are multiple ways of managing the departure queue area, namely FIFO queues where 
aircraft need to be assigned to the queues, FIFO queues where queues are pre-assigned, and no queues. The 
formulation presented here addresses all these configurations. Thus, runway usage is modeled as a double 
assignment problem; aircraft expected to depart within the planning horizon are assigned to a departure queue and 
then all aircraft, either arrivals or departures, are assigned a runway usage time based on various separation 
constraints. The approach is generalized for use in a variety of situations and allows for aircraft prioritization based 
on operational considerations. 

Besides the above three cases for modeling, scheduling can be directed towards three basic objectives (or any 
combination thereof): throughput, system delay and maximum delay. Throughput is linked to the capacity of the 
runway and for scheduling purposes can be defined as the assigned take-off or crossing time for the last aircraft 
being scheduled. System delay is linked to efficiency and is the total time spent by all aircraft while waiting to take- 
off or cross the runway. Maximum delay is linked to equity and is the maximum waiting time for any aircraft while 
waiting to take-off or cross. The MILP described in this paper can address each of these three objectives as 
described in the next section. 
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Figure 1 : Example departure queue structure and arrival crossings for runway 17R at DFW 


1 1 1. Mathematical Details 


A. Basic Mathematical Formulation 

Consider a set of flights F , which represents the aircraft expected to use the departure runway (either take-off or 
crossing) within the planning horizon. Let A and D be two subsets of F, denoting arrival flights and departing flights 
respectively. Each member of the set D will have an aircraft type and destination departure fix, which would govern 
the separation from preceding and succeeding aircraft. Each member of A will have an aircraft type and the assigned 
crossing. Further, each aircraft in F will have a time of earliest runway usage based on when it enters the queuing 
area (un-impeded, or assuming no other aircraft is present) or is available to cross. Let this time be df for flight 
/ (/ E F), and let the set F be sorted in ascending order of df. In other words, if i,j E F then i > j implies a t > dj. 
Each aircraft may have a window of runway usage, and let us denote this with Af. This time window could be strict 
(that the aircraft has to depart before df + Af) in which case smaller values of A j could lead to infeasibility, or the 
window could be lenient, with a penalty for violation. In the formulation presented here the time window is modeled 
as a strict constraint. With the small modification, this formulation can be changed to accommodate lenient, 
penalized time-windows too. These windows can be used as a means to prioritize some aircraft. Of course, 
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prioritization is optional and the windows could be sufficiently large to avoid any prioritization. Let P denote the set 
of positions these aircraft can occupy in the “usage sequence” from the runway. Thus, \F\ = \P\, or the number of 
such positions in the output sequence are equal to the input number of aircraft. The rest of this section describes the 
notation, summarizing the above set of parameters and describing the variables. This is followed by the complete 
mathematical formulation of the MILP and an explanation of all the expressions. 

Parameters 

F is the set of all incoming flights, sorted in ascending order of earliest un-impeded time at runway 

A is the set of all arrival flights, A c F 

D is the set of all departure flights, Dcf 

P is the set of positions that flights can take in the departure sequence on the runway, \F\ = |P| 

Q is the set of all departure queues that the departure aircraft can be assigned to 

C is the set of all crossings for the runway, and each arrival flight is on one of these crossings 

df is the earliest runway usage time for flight / (/ E F) 

A f is the time window of runway usage for flight / (/ E F) 

hf c is the additional waiting time before arrival flight / crosses due to the usage of crossing c, (/ E A,c E C) 

Wf C is a binary parameter, and is 1 if c is the crossing for arrival flight /, zero otherwise (/ E A,c E C) 

dij is the minimum separation in runway usage times of flight i and j, if flight i follows flight j ( i,j E F) on 
the runway. This depends on a variety of factors, as described below. 

The parameter d^ defines the required separation between the two aircraft, and depends on whether the aircraft are 
arrivals or departures, and in the case of arrivals, the assigned arrival crossing. Let 

e\j be the required separation in runway usage times of departure flights i and j, if flight i follows flight 
j ( i,j E D ) on the runway. 

efj be the required separation in runway usage times of departure flight i and arrival flight j, if flight i follows 

flight j ( i E D,j E A ) on the runway 

efj be the required separation in runway usage times of arrival flight i and departure flight j, if flight i follows 

flight j ( i E A, j E D ) on the runway 

g c be the required separation between two arrival flights when both are assigned to crossing c, (c E C) 


Note that e\j is the minimum separation between the two departure aircraft considering all the separation criteria 
including wake-vortex separation, miles-in-trail restriction etc. Further, depending on the assigned crossing, there 
could be additional delay ( hf C )• Incorporating all these cases, the parameter d t j is defined as below: 



(Uj e D) 

(i E D,j E A ) 

O' e A, j E D) 

0 i’j eA ’ w i,c = Wj,c = l,c EC) 
0 ij e A ' w i,c + w j,c < 1 ,c EC) 


( 1 ) 


Variables 
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is a binary variable, and is 1 if flight / (/ G F) occupies position p (p EF), zero otherwise 

is a binary variable, and is 1 if departure flight / (/ED) is assigned to queue q (q E Q), zero otherwise 

is the runway usage time of aircraft assigned to position p 


Vf.P 

x f,q 


Mathematical Formulation 


Minimize System Delay: minimize 


Y tp Y Uf 

pep feF 


Maximize Throughput: minimize t\ P 


Minimize Maximum Delay: minimize max 

pep 



( 2 ) 

(3) 

(4) 


such that 

Y yfp = 1 V f eF (5) 

pep 

Y y f’ p = 1 VpeP ( 6 ) 

feF 

Y X f« = 1 VfeD ( 7 ) 

Q E Q 

VpeP (8) 

feF 

t v <Y y f , P ( a f + b f ) VpeP 

(9) 

feF 

t p > t v < + (y j p + y j v ' - 1 )d u Vi,j e F, p,p' e P,p > p' (10) 

^(y;,u)w-^(yj>)v > \P\{x Uq +x jiQ - 2) VqEQ, ij ED,i>j (11) 

ueP vep 

^(yu) u ^ w t)C =Wy )C = l, cec, i > j, i,j e D (12) 

ueP i?ep 


As stated before, each departure aircraft is assigned to a departure queue and then assigned a departure time. The 
time assignment is done by first assigning each aircraft a position in the output departure sequenc e (if it is a 
departure aircraft), and then maintaining the separation between consecutive positions. The expressions | (2) | to [~( 1 1 ) | 
can be interprete d a s fol lows: 

Eq. (2)]|(3) and |(4)| are the three objective functions labeled according to the pertinent objective. Eq. |(2)| is the 
system delay objective, and minimizes the total time spent by each aircraft in departure queue (or crossing for 
arrivals) beyond the earliest, unimpeded runw ay usage time. Eq. |(3)| is the throughput objective that minimizes the 
take-off time of the last position, t\ P \. Eq. (4) is the maximum delay objective, which minimizes the maximum time 
spent by any aircra ft in the queue area. 

Constraint ! (5) | mandates that each aircraft will occupy just one position in the output sequence. 
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Constraint ! (6) | is similar to |(5)J and represents the fact that each position in the output sequence must be occupied 
by just one ai rcraft . 

Constraint ) (7) |e nsu res th at each departure aircraft can occupy just one departure queue. 

Constraints ! (8) | and | (9) | ensure that if aircraft / is assigned to position p, then the runway usage time for position 
p is within th e time window of aircraft /. 

Constraint ! (10) | ensures that the runway usage times for any two positions are sufficiently spaced out to ensure 
the separation criteria flights assigned to the positions. If flight i is assigned to position p (y i p = 1) and flight j is 
assigned to position p — k {y JiP -k — l), then t p > t p _ k + dy, otherwise the constraint is “inactive,” or dominated 
by the constra int for the pair of aircraft that are assigned to these positions. 

Constraint ! (11) | ensures the FIFO structure of each departure queue. The two terms on the left hand side of the 
inequality give the difference in the position of aircraft i and (i > j). If both the aircraft are assigned to the same 
queue {x iq = Xj q = l), then the difference in the positions should be at least one, with i occupying a position later 
than j since i entered the queue later than j. Similarly, constraint ( 1 2) ensures the FIFO structure of the arrival 
crossings. 

To address scenarios with pre-determined departure queues, additional constraints can be included in the 
problem, such as Xf q = 1 for the relevant queue and departure flight. The rest of the formulation remains the same. 
In fact, for these scenarios, Xf q is no longer a variable but an input parameter, and thus the problem size is reduced. 

For scenarios with no FIFO queues, essentially there is no d epart ure queuing structure. This can be modeled by 
eliminating the queuing variables Xf q and related constraints (7) and (11) from the formulation. The reduced 
problem will have fewer variables and constraints, and depending on the problem characteristics, it is possible that 
the reduced problem is computationally more efficient. 

As stated before, the above formulation models the time window as a strict constraint. This formulation can be 
changed to accommodate “soft” time-windows by defining an additional linear variable s p and adding constraints 
s p > t p — Y*fyf, P (af + A/) and s p > 0. Then s p can be multiplied by an appropriate cost of violating time windows 
and added to the objective function. 

B. Modifications for Computational Improvement 

Computationally, the formulation presented above does not pe rform very well. One of the reasons for poor 
performance is the number of constraints presented in constraint ! (10) J For a 25 aircraf t prob lem, this results in 
almost 165,600 constraints with binary and linear variables. For this reason, constraint ( 1 0) is re-formulated as 
below: 

i*j 

t v - t v < > Y y Up d Uj - (l - y JiPI ) (max d Kl - min d w ) Vp,p' E P, j E F (13) 

iEF 

The above constraint is an aggregated form of constraint (10) and reduces the large number of separation 
constraints while retaining the same bounds. For a 25 aircraft problem, th is redu ces the number of constraints to 
6,900. This reformulation is similar to the one presented in Ref. 9. Constraint ) (1 l)| can be reformulated as : 

VEP UEP 

*i,q + *j,q + Y yiv + Y y i’ u ~ 3 VU)£D,i>j,\/pEP,VqEQ (14) 

v=p..\P\ u=l..p 

The above constraint states that if both departure aircraft i and j are assigned to the same queue q and a t > a ; -, then 
for any position p either i takes the position p or greater than p, or j takes the position p or a position less than p. 
Since the abo ve con straint does not incl ude no n-unitary coefficients for the binary variables, it gives better bounds 
than constraint ) (1 1)J Similarly, constraint ! (12) | can be reformulated as: 

VEP UEP 

Y yi ’ v + Y y ’’ u - 1 v i’j e A ’ 1 > j' v V 6 p - w i,c = w j,c = 1 ,c EC (15) 

v=p..\P\ u=l..p 

Other constraints presented in Ref. 9 can also be used here, with varying degrees of effect. However, the 
constraint based on the best solution for all = 0 is not apparent because of separation being governed by crossing 
assignment of arrival aircraft. Generating such a “best solution” for better lower bounds is a direction for future 
research. 
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IV. Results 


A. Benefits over F i rst-Come- F i rst-Se rve 

The results from the runway scheduling using the above MILP are compared with a basic first-come-first-serve 
(FCFS) rule to e xami ne the benefits of the scheme. The benefits are presented for two of the three objectives listed 
earlier in Section flllj throughput and system delay. As discussed in Ref. 9, system delay and throughput seem the 
primary objectives used by the controllers today and using the maximum delay objective results in large variations 
in system delay and throughput values. The effect of using these two objectives on maximum delay is investigated 
later in this section. 

For the purpose of comparison, randomly generated problems comprising 15 departure aircraft and 10 arrival 
aircraft were used. The description of the problems is given below in the next paragraph. The MILP was solved 
using IBM ILOG CPLEX Optimizer 4 , a commercially available optimization software package, and the MILP was 
solved for an optimality gap (% gap in lower and upper bound) of 1.0%. Again, of the three possible configurations 
of the departure queue area as described in Ref. 9, results are provided only for the case with FIFO queues, where 
the queue assignment is not pre-determined. Based on operations at DFW, three queues and four arrival crossings 
were used in the analysis. 

To generate the problems, data was collected using the Surface Operations Data Analysis and Adaptation 
(SODAA) tool 5 * * . This tool analyzes the Surface Management System 11 (SMS) generated log files, which contain data 
from multiple sources, including air carriers, the Enhanced Traffic Management System (ETMS), and Airport 
Surface Detection Equipment, Model X (ASDE-X). Four aircraft weight classes were used for which FAA mandates 
minimum separation criteria based on wake vortex separation 12 . This separation is given in terms of distance, and 
was converted into time -based separation using average runway occupancy times and average horizontal and ascent 
speed for each airc raft weigh t class based on actual surface data at DFW airport. The resulting time -based separation 
matrix is given in | Table l] Not e that the values for arrival aircraft are dependent on the assigned crossing, as 
discussed before in equation [(T)] Thus, the values need to be augmented with 0. It was assumed that 0 = 0 for the 
crossing closest to the departure queue, and 0 = 3, 6 and 9 seconds for the remaining three crossings, with the value 
increasing with increasing distance from the departure queue. For two successive arrival aircraft from the same 
crossing, a required separation of 40 seconds was added. This is due to the practice of controllers providing 
sufficient spacing between crossings from the same queue to avoid arrivals getting stuck on the runway due to 
congested taxiways immediately after crossing. 

The earliest runway usage times (a^) were uniformly distributed, and the aircraft type is randomly assigned 
using a uniform distribution (resulting in approximately 25% aircraft of each weight class). Also, the A^ values were 
set to infinity to allow for maximum re-sequencing. Setting a lower value would be equivalent to providing upper 
bounds to the maximum delay objective since A ^ is the maximum allowed delay for flight /. 


Table 1 : Aircraft Types and Minimum Departure Separation (in seconds) 



| Figure 2| shows the improvement in objectives over FCFS for 50 randomized problem sets with uniform aircraft 
mix (i.e., 25% each of small, large, heavy and B-757). The different lines in the figure represent the percentage 
improvement over FCFS for the two objectives. Both the percentage gains and the absolute benefit of all cases for 
both objectives are presented. 


4 http://www.cplex.com 

5 http://www.mosaicatm.com/Projects/SurfaceOperationsDataAnalysisandAdaptation/ 

7 

American Institute of Aeronautics and Astronautics 





The gains in throughput are relatively small as compared to the potential reduction in system delay (average 
throughput improvement is appr oximately 9%, whereas average system delay reduction is approximately 50%). 
However, the potential benefits in Figure 2| are for a uniformly distributed aircraft mix (i.e., approximately 25% of 
each of the four classes defined in Table 1). Although it can be argued that the potential aircraft mix in the future 
would change, benefits on current aircraft mix should also be explored. For this purpose, the data analysis from 
SODAA was used to determine an approximate mix of aircraft at DFW: 2% small, 88% large, 5% heavy and 5% B- 
757. In Figure 3, the potential benefits for the current aircraft mix are provided. Again, these results are for 50 
problems of 25 aircraft with 15 departures and 10 arrivals. The observed gains in throughput (average improvement 
of approximately 7%) as well as system delay (average improvement of approximately 44%) were slightly less than 
those in uniform aircraft mix. Again, the gains in system delay were significantly more than throughput. 


System-delay improvement over FCFS 
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Figure 2. System delay and throughput improvement over FCFS for uniform aircraft mix 
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Figure 3. System delay and throughput improvement over FCFS for current aircraft mix 


B. Objective Comparison 

In the previous section, potential benefits from scheduling directed towards an objective in isolation were 
presented. However, it is necessary to understand the relative effect of using one objective on the other metrics. 
Scheduling using just one metric as objective can have adverse effects on other metrics. For example, if scheduling 
is done based on the maximizing throughput, the optimal solution could result in five aircraft being made to wait to 
let another aircraft take-off first for slightly better throughput, resulting in large delays and reduced efficiency. In 
this section, the relative merit of choosing system delay versus throughput is first compared, and then the choice of 
maximum delay as an objective is discussed. These comparisons are done for both uniform and current aircraft mix. 
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System Delay vs Throughput 
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Figure 4. Comparing system delay and throughput objectives 


The 100 problem sets (50 for uniform mix and 50 for current mix) were used, and | Figure 4| shows the result of 
this comparison. Each dot on the figure represents one problem set. The horizontal axis represents the percentage 
deviation from optimal throughput while minimizing system delay and the vertical axis represents the percentage 
deviation from least system delay while optimizing throughput. The figures show very small deviation in throughput 
when system delay is optimized (at most around 6%), but a very large deviation in system delay for optimal 
throughput. This indicates that multiple throughput-optimal solutions are present, or multiple throughput 
maximizing solutions exist in the close vicinity of the optimal solution. Searching for the least system delay solution 
in this neighborhood or using system delay as the objective itself might be more relevant from the point of view of 
future implementation. 

Ref. 9 shows that system delay is more sensitive to the choice of objective function (a statement strengthened by 
Figure 4j . Results from using maximum delay as the objective function showed large deviations even for the 


formulation presented here. Maximum delay becomes important for equity reasons: throughput or system delay 
improvements at the cost of significant delays to certain aircraft might not be an acceptable solution. To see the 
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effect of using throughput and system delay (as objective functions) on maximum delay, we compared the maximum 
delay from the optimal solutions for these objectives with the maximum delay from the FCFS solution. It can be 
argued that the FCFS solution would provide an “upper bound” to the de lay of each ai rcraft, with no aircraft 
experiencing more delay than what it would face in the FCFS solution. | Figure 5 and Figure 6 show these 
comparisons for uniform mix and current mix respectively. The figures show large variation in maximum delay 
between system delay or throughput optimal solution and the FCFS solution. To conclude, even though system 
delay yields an almost optimal solution for throughput, using it as an objective in isolation may not result in 
acceptable solutions. One possible method would be to use system delay as the objective with the FCFS maximum 
delay for each flight being used as a constraint through Af. 




Figure 5. Comparing max delay in FCFS solution with optimal system delay and optimal throughput solution 

for uniform aircraft mix 




Figure 6. Comparing max delay in FCFS solution with optimal system delay and optimal throughput solution 

for current aircraft mix 


C. Computational Performance 

The abov e for mulation does not perform very well computationally, even while using the improvements presented 
in section fin] The average solution time for each run above was more than 30 minutes. One of the reasons for the 
long computation time was the lack of good lower bounds on t p , demonstrated by the fact that the solution at the 
end of 2 minutes computation time was almost always the optimal solution. Enhancing the computational 
performance of the above method is a direction for future research. 
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V. Conclusions 

This paper presented a mixed integer linear program for deterministic runway usage scheduling. The model is 
generic and can be used for a variety of cases with different methods of handling the queuing area. The MILP 
explicitly considers separation criteria along with additional constraints and includes an optional prioritization 
scheme for relevant aircraft. Multiple objectives can be used and simulations indicate substantial benefits over a 
basic first-come-first-serve rule. Results show the largest benefits in reducing system delay, with system delay being 
the more sensitive objective compared to throughput. Minimizing system delay alone also provides benefits in 
throughput. Computational improvements to the basic MILP are also provided; however, in almost all cases the 
large solution times are due to poor bounds, i.e. the optimal solution was found fairly quickly and a lot of time was 
spent in proving optimality for this solution by changing the lower bounds. 

Since the model solves the deterministic problem, the logical next step is to address uncertainty in the earliest 
runway usage time while solving for larger time periods. One way to do this would be to use rolling planning 
horizon methods, and the use of the above model in such a scheme needs to be studied. Further, improvement in 
computational times is ongoing research, to make the model more relevant for implementation in a decision support 
tool for controllers. 
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